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Abstract. We present a simple approach for finding A-soliton solution and the 
jj^ , corresponding Jost solutions of the derivative nonlinear Scrodinger equation with 

P I ' nonvanishing boundary conditions. Soliton perturbation theory based on the inverse 

(-h . scattering transform method is developed. As an application of the present theory we 

consider the action of the diffusive- type perturbation on a single bright / dark soliton. 



PACS numbers: 05.45.Yv, 52.35.Bj, 42.81.Dp 
1. Introduction 



The derivative nonlinear Schrodinger equation (DNLSE) (a = ±1) 

id t u + d 2 u + iad x (\u\ 2 u) = (1) 



has many physical applications, and, probably, the most important are in plasma 
physics and in nonlinear optic. First, equation (pQ) describes modulated small-amplitude 
nonlinear Alfven waves in a low-/5 (the ratio of kinetic to magnetic pressure) plasma, 
propagating parallel [H El [3] or at a small angle jU El [7] to the ambient magnetic 
field. The DNLS equation also describes large-amplitude magnetohydrodynamic waves 
in a high-/? plasma, propagating at an arbitrary angle to the ambient magnetic field 
[8]. In these cases u denotes the transverse magnetic field perturbation normalized by 
the ambient magnetic field, where t and x are normalized time and space coordinates, 
respectively. Second, the DNLSE is related to the modified nonlinear Schrodinger 
equation (MNLSE) 



id T *p + G -d\^ + is0f(|V>| V) + lVfV> = (2) 
by a simple gaugelike transformation [9] 

V>(£,t) = u(x,t)exp j^ji + ^xj , (3) 

where a = ±1 corresponds to the abnormal (normal) group velocity dispersion (GVD) 
region, f = fax) /(2s) + t/(2s 3 ), r = fat)/(2s 2 ). In turn, the MNLSE describes the 
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propagation of ultrashort femtosecond nonlinear pulses in optical fibers, when the 
spectral width of the pulses becomes comparable with the carrier frequency, and, in 
addition to the usual Kerr nonlinearity, the effect of self-steepening of the pulse should 
be taken into account. In this case, u is the normalized slowly varying amplitude of the 
complex field envelope, t is the normalized propagation distance along the fiber, and x 
is the normalized time measured in a frame of reference moving with the pulse at the 
group velocity 

Equation (TjQ) is completed by the boundary conditions: vanishing (u — > as 
|x| — > oo) or nonvanishing (\u\ — > p = const as |x| — > oo) at infinity. In both cases 
the DNLSE is integrable by the inverse scattering transform (1ST) [TUJ [HI [12J, [13] , and 
admits iV-soliton solutions [Tij . 

The nonvanishing boundary conditions (NVBC) are important in physical 
applications. For example, in space plasma physics the vanishing boundary conditions 
(VBC) are relevant only for the case of propagation of Alfven waves strictly parallel to 
the ambient magnetic field. In nonlinear optics the NVBC can support propagation of 
dark solitons in both normal and abnormal GVD regions [15J. Unlike the nonlinear 
Schrodinger equation or the DNLSE with VBC, the DNLSE with NVBC admits 
simultaneous generation of breathers (solitons with internal oscillations) and one- 
parametric (nonoscillating) bright and/or dark solitons [16]. 

The 1ST formalism for the DNLSE with NVBC is much more complicated from the 
one for VBC. Analytical properties of the Jost solutions in this case are formulated on the 
Riemann sheets of the spectral parameter [TT] and the corresponding direct and inverse 
scattering problems are rather involved. Recently, Chen and Lam [15] developed the 1ST 
for the DNLSE with NVBC by introducing an affine parameter to avoid constructing the 
Riemann sheets. Both approaches, however, encounter a difficulty when finding exact 
explicit iV-soliton solutions. The reason is that the resulting solution u containes the 
phase factor exp(ir] + ), where i] + is some definite integral from \u\ 2 . Thus, the solution 
is written in an implicit form and only modulus of the solution can be obtained in that 
way. Though for simple one-parametric soliton solutions the phase r] + can be calculated 
by direct integration, this procedure is obviously impracticable for iV-soliton solutions. 
Instead, tricks leading to the explicit expression for r] + were used in some particular 
cases: for the two-parametric one-soliton breather solution [T5] . and for the iV-soliton 
with purely imaginary discrete spectral parameters (i. e. for pure bright and/or dark 
solitons) [TTl [T8]. Another approach based on Darboux/Baclund transformations was 
developed by Steudel [H]. Apparently, Steudel was the first to obtain exact iV-soliton 
solutions with explicitly calculated phases for the DNLSE with NVBC. 

From the practical point of view, the completely integrable DNLSE ([1]) is an 
idealized model. In many physical applications, additional terms are often present in 
the DNLSE. They can include effects of the third-order linear dispersion, dissipation, 
influence of external forces, etc. . These terms violate the integrability, but being small 
in many important practical cases, they can be taken into account by perturbation 
theory. The most powerful perturbative technique, which fully uses the natural 



N-soliton solutions and perturbation theory for the derivative nonlinear S eroding er equation with nonvanish; 



separation of the discrete and continuous (i.e., solitonic and radiative) degrees of freedom 
of the unperturbed DNLSE, is based on the 1ST. While the IST-based perturbation 
theory for the DNLSE with VBC was developed long ago [T9j , the analogous theory for 
nonvanishing boundary conditions was absent. 

The aim of this paper is twofold. First, we present a relatively simple approach 
for finding exact explicit (i. e. with the phase) N + M-soliton (N breathers and 
M "usual" bright and/or dark solitons in asymptotics) solutions of the DNLSE with 
NVBC and show that these solutions can be obtained without determining the phase 
factor exp(ii] + ). Thus, exact exotic solutions, describing, for instance, collisions between 
breathers, as well as collisions between pure bright/dark solitons and breathers can be 
written. Simultaneously, unlike the purely algebraic approach jUj based on Darboux 
transformation, the corresponding Jost solutions can also be obtained. A second aim 
is to present perturbation theory for solitons of the DNLSE with NVBC. We derive 
evolution equations for the scattering data (both solitonic and continuous) in the 
presence of perturbations. As an application of the present theory we consider the 
action of the diffusive-type perturbation on a single bright / dark soliton. 

Without loss of generality, we will consider the NVBC in the form 



We also put a = 1, since the case a = — 1 can be obtained from the former by a 
transformation x — > —x. 

The paper is organized as follows. In section 2, we review a theory of the scattering 
transform for the DNLSE with NVBC. In section 3, we present the procedure to 
construct the general explicit iV + M-soliton solution. Integrals of motion are obtained 
in section 4. The perturbation theory and its application are considered in sections 5 
and 6, respectively. The conclusion is made in section 7. 

2. Inverse scattering transform for the DNLSE with NVBC 

In this section we present the theory of the scattering transform for the DNLSE with 
NVBC, following [15] with some modifications and specifications. Equation ([ID can be 
written as the compatibility condition 



u — > pexp(±2z6 l ), 



as 



x 



±oo. 



(4) 



d t U-d x V+[U,V} = 0, 

of two linear matrix equations (Kaup- Newell pair) [10] : 

d x M(x,t, A) = UM(x,t, A), 
d t M(x, t, A) = VM(x,t, A), 



(5) 



(6) 
(7) 



where A is a spectral parameter, and 






V = -2i\ 4 a 3 + 2A 3 Q - i\ 2 Q 2 a 3 + AQ 3 - i\Q z a 3 



(9) 
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Consider the linear problem for some fixed t. In terms of the matrix U boundary 
conditions (j4]) can be rewritten as lim U(x, A) = ^(A), where 

X— >rboc 

u±=( - l>? " Xe± ™\ (10) 

Asymptotic solutions of ([6]) E ± satisfy 

d r E ± = U ± E ± . (11) 



The double-valued function K(X) = Aa/A 2 + p 2 appears in the matrices E ± } and the 
analytical properties of solutions of Eq. ([6]) are formulated on the Riemann surface 
determined by the function K(X). The Riemann surface S consists of two sheets S + 
and S~ of the complex A-plane with branch cuts on the image axis from — oo to — ip 
and from ip to oo. It is convenient to introduce an affine parameter ( by a change of 
variable [15] 

A(C) = i(C-£V (12) 



c 

This transformation maps the sheets S ± onto Im£ > and lm( < respectively and 
the real axis on the complex A-plane into the real axis on the £-plane. Under this, 

E ± (x,() = ^ J p/( . "^ /C \ e- lk ^ x (13) 

where the single-valued function k(() is determined by 

hq = 1(c+£)ko (14) 



2 v C 

For lm( = denote by M ± (x,C) the 2x2 matrix Jost solutions of ([6]), satisfying the 
boundary conditions 

M ± -> E ± (x, C), as x -> ±oo. (15) 

The corresponding integral equation can be obtained from ([6]) and (TT5T) 



/•±oo 

r±l • x 5 - E ± (x, A) =F iA / E ± (x - y, X)Q(y)M ± (y, A) dy. (16) 

J X 

The matrix Jost solutions M ± (x,C) can be represented in the form M~ = ((p,0) and 
M + = (-0,-0), where 9? and ip are independent vector columns. The monodromy matrix 
S(() relates the two fundamental solutions M~ and M + : 

M-(x,() = M + (x,()S((). (17) 
The Jost coefficients are defined by 

ip = aip + bip, (18) 

(p= -aip + bip, (19) 



M ± (x,A) 
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so that the monodromy matrix is 

*> 

where ad + bb = 1. Matrices and S have the parity symmetry properties 

S(C) = CT 3 S(-C)a 3 , M ± (C) = a 3 M ± (-C)a 3 , (21) 

and the conjugation symmetry properties 

S(C) = <t 2 S*(C)o- 2 , M ± (C) = a 2 M ± *(C)a 2 , (22) 

where a 2 and 03 are Pauli matrices, so that \a\ 2 + \b\ 2 = 1. In addition, since the 
scattering problem ([6]) possesses symmetry with respect to the inversion £ — > p 2 /(, the 
important involution properties are valid: 

M ± (x,p 2 /C) = (C/p)^ 3 M ± (x,C)a 2 , (23) 
S(p 2 /C)=a 2 S(C)a 2 = S*(C) (24) 
It follows from (fTTl) that 

a(0 = A- 1 (C)det(^V), (25) 
where we have introduced the notation 

A(C) = detM ± = l + p 2 /C 2 - (26) 

Columns ip(x,() an d i/j(x,() turn out to be analytically continuable to lmk(() > 
(i. e. to the first and the third quadrants of the complex £-plane), while (p and 
•0 are analytically continuable to lmk(() < (i. e. to the second and the fourth 
quadrants) [TT1[T5]. Then, the coefficient a(Q is analytically continuable to lmk(() > 0. 
The analytic function a(£) may have zeros (i,...,(n m the region of its analyticity 
Im fc(C) > 0. Equation (|25l) then shows that the columns i/j and <p are linearly dependent 
and there exist complex numbers bi, . . . , fejv such that 

ip(x,Q = bjip(xXj), (27) 

and, similarly 

<p(x,Q) = -b* j $(x,q). (28) 
The standard analysis of (1TB]) yields the asymptotics at |£| — >■ oo 

IK*, C)e- ifc( ^' - (~™ /C ) e*^ + 0(1/|C| 2 ), (29) 



where 



<p(x,()e ik{C)x - (_-^ /c )^"- 0) + 0(1/|C| 2 ), (30) 
T) ± = ±~ / (p 2 - M 2 ) dx. (31) 



2 
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(p(x, ()e ik{0x - 



As \C\ — > 0, we have 

-ip/C 
u*/p 
u/p 

-ip/C 

It then follows from ( 1251) that asymptotics of a(Q are 



-i(j7+-6 



-0(1), 
0(1), 



a (C) — y exp(ii] — 2i9), as \(\ — 
a(() ^> exp(-ir] + 2i6), as |C| 



oo 



where 



1] = 1] + 1] 



1 f°° 

2 y - n 2 )^- 



(32) 
(33) 

(34) 
(35) 

(36) 



Zeros of a(£) in the region of its analiticity (i. e. to the first and the third quadrants 
of the complex £-plane) are not independent due to the symmetry properties (12T1) . (1221) 
and (1241) [T5] . If Q is a simple zero of a(() in the first quadrant, outside the p circle, 
then — Q (in the third quadrant), p 2 /Q (in the first quadrant and inside the p circle) 
and —p 2 /Q (in the third quadrant and inside the p circle) are also simple zeros of a(£). 
There are only two zeros for each j if Q lies on the p circle: (j and —Q. Thus, one can 
consider zeros Q lying only in the first quadrant outside and/or on the p circle. These 
zeros can be parametrized as Q = pexp(7j + i/3j), where jj ^ and < f3j < tt/2. 
In what follows we assume that in the first quadrant M zeros lie on the p circle and 
N zeros lie outside the p circle so that j = 1 . . . M + N. Using the asymptotics ( 1341) . 
(135]) and standard methods of the Hilbert transform theory [20] in conjunction with the 
properties (l2~Ij) . ( |22|) and ( l24|) . one can express the coefficient a(() in terms of its zeros Q 
in the first quadrant outside and/or the p circle, and the values of |6(C)| on the contour 
T = [0, -oo] U [0, oo] U [zoo, 0] U [-zoo, 0] 



«(C) 



x exp 



1 

2tt7 



(C 2 - C 2 ) (C 2 



p7c; 2 ; 



N 

n<c 2 -c* 2 ) (c 2 -pvci) 



M f/-2 A2 



i« 5 



G 2 ! 



Setting ( 



P-C 

in (1371) and comparing with (132|) . we get 

M N 

7 1 = 26-2j2^g(k - 4^ argC, + — 

fe=l j=l 



HHp)? 
p 



-dp,. 



(37) 



(3? 



The potential in the general case is 

f N 

- 2pe-'^ + -e) 

[j=i 



u(x) = pe 



-2i(r]+-e 



M 



1 



(39) 
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where Cj = bj/a'- with a'- = da/d^\^ = ^. For the compatibilty with the second Lax 
equation ([7]), the Jost solutions obtained from (jfJJ) should be multiplied by a t-dependent 
factor h((,t) = exp[-zfi(C)*], where = [2A 2 (C) - p 2 }k(() [15J: 

^(x, C, t) = h(C, t)$(x, C), ^(x, C, t) = h-\c, t)^(x, C), (40) 

f(x,C,t) = h(t,t)<p{x,t), <p(x,t,t) = h-\C,t)<P(x,0- (41) 
Dynamics of the scattering data turns out to be trivial 

a(C,t)=0, (42) 

6(C,t) = 6(C,0)exp[2ifi(C)t], (43) 

6 i (*) = 6 i (0)exp[2ifi(C i )t],. (44) 

3. The Jost solutions and the potential in the reflectionless case 

An important particular case is that of the reflectionless (solitonic) potentials u(x) when 
b(t, C) = as a function of ( for some fixed t. It then follows from (l3"Tj) and fl38l) that 

ofn n C(c 2 -ci) A c; 2 (c 2 -c|)(c 2 -p7c; 2 ) 

M ^(C 2 - C 2 ) }J C|(C 2 - Cf) (C 2 - pVC 2 ) ' 1 j 

which extends to Im( 2 < as a meromorphic function. One also sees that a(t,() = 
l/a(t,()- Since S(t,() is diagonal in this case, it can be factorized in such a way 
S ~(£) = S + (()S((), that the Jost solution matrices M ± is expressed through a common 
solution matrix A(x, () 

M ± (x,()=A(x,()S ± ((), (46) 



where 



with 



5 " = n (c^c!) n c;(c 2 - c 2 )(c 2 - p 4 /c; 2 ) ' (48) 



and 



S- = a 1 S + a 1 , (49) 

where ci is Pauli matrix. One can see from (j4Tj) and f|49l) that S + (p) = S~(p) if M is 
even, and 5* + (p) = S~(p) if M is odd. Therefore, from Q4"5T) we get 

Af+(z,p) = (-l) M M"(x,p). (50) 

On the other hand, since A = corresponds to £ — ±p, we have from ( JTBl) and ( 1T61) 



.U ± (,-./,) = ^ ± f.r. = ( ^ ^ J . (51) 
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It then immediately follows from (joT)j) and (joTj) that 9 = ± irn if M is even, and 
# = ±7r/2±7rn if M is odd (n is an integer). Thus, we established the following important 
fact: the total phase shift 46 in the iV-soliton solution is zero (or an integer times 27r). 
Note, that authors of [5J [15] showed that 9 = for the particular case N = 1 considering 
an explicit, rather complicated expression for the one-soliton breather solution. From 
f l46|) - fj48|) and (l5Tj) one can also obtain 

a (x , „j = p-(- ffl) M ( \ ) m - w n '^^i?/" 10 ' 4 ' - 

V / fe=l j=l l^'l 

As follows from (I2"3l) . ([27)1 and (|28|) the columns of A(x, £) satisfy the relations 

A i{x, Cj) = bjA 2 (x, Q, (53a) 
A 2 (x,C;) = -b* j A 1 (x,C]), (536) 

A 1 (x, P 2 /c;) = ^ 2 (x,p 2 /c;), (53c) 

A 2 (x,p 2 /Q) = -Mi(^>P 2 /Ci), (53d) 

for all j = 1, . . . , M + N. For M zeros lying on the p circle, equations (l53al) - (153aT) 
become 

A 1 (x,Q) = b j A 2 (x,Q), (54) 

A 2 (x,(*) = -b 3 A 1 (x,q), (55) 

where the coefficients bj are real. One can see from ([32]) and ()4"6l) that A(() is analytical 
in the whole ( plane, except for the point ( = 0, where off-diagonal elements of 
A(() exp[ik(()xos} have simple pole. Thus, the matrix C^(C) ex p[^(C) 2;cr 3] is analytical 
in the whole ( plane. It then follows from (|29|) and (146]) that diagonal and off-diagonal 
elements of the matrix C^(C) ex P[^(C) a;<J 3] are polynomials in ( of degrees 4iV + 2M+ 1 
and 4N + 2M, respectively. In addition, from ([2~T|) and (piET) one sees that the diagonal 
elements of A are even in £, while the off-diagonal ones are odd. This means that we 
can write 

A(x,()e^=[ 21-, ^!(o) )+YA 2p ~ 1 \ C fl A l ), (56) 




22 



where L = 2N + M, and A™ are some unknown functions of x. Setting ( = p in (1561) 
and comparing with 




4? 4V 1 

.(0) 1 ,(0) }^A^I J ' I jto _4(P) 
^21 P A 22 / „=! \ A 2l P A 22 



P 2 "(-<71) M 



n ^g±f^M . (57) 



The unknown coefficients Am^i, t) with p = 0,...,L are determined uniquely from 
([3]) and (1571) . Indeed, the first row of ([3]) and (l57j) is a linear algebraic system of 

g and 4^ 



2L + 2 equations in 2L + 2 unknowns, the coefficients 4^ and with p = 0, . . . , L 
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Figure 1. (a) Bright and dark soliton amplitudes versus the parameter f3 for p = 1. 
The dotted line and the solid line up to the bifurcation point correspond to the bright 
soliton. The solid line corresponds to the dark soliton. The bifurcation point is at 
Per = 7r/6 and |ito| = P- (b) Dark soliton profiles for different f3 and p = 1. Curves 1 
("black" soliton), 2 and 3 correspond to [3 cr = ir/6, (3 — 0.3, and [3 = 1.1, respectively. 



Likewise, the second row of ([3]) and fl57|) is the system for determining and 
with p — 0, . . . , L. By direct substitution one can check that ( |56i) is compatible with 
© and dUD if and only if 

i \x, t) 

This formula reconstructs u(x, t) from the discrete scattering data {(j{t)}, {bj(t)} in the 
case when b(t, () = and it gives N + M-soliton solution of ([1]). An explicit form of the 
solution can be easily written in terms of the determinants of corresponding matrices. 
Equations ()46l) and ( 1561) determine the N + M-soliton Jost solutions. 

As the first example, let us consider the simplest case when the function a(() has 
one simple zero (i in the first quadrant of the complex (-plane on the p circle (i. e. 
M = 1, N = 0) so that d = pexp{ij3i) with < fix < tt/2. Taking into account Eq. 
()44l) . we have b\ exp(2i/cix) = eexp(— z), where z = k (x — vt — x$) with 

k = p 2 sin(2/?i), v = 2p 2 - p 2 cos(2/?i), (59) 

and without loss of generality one can set e = ±1. Determining A^J and AJq from ([5 



and solving a system of two linear algebraic equations for A^ 2 and A^ from ( 1541) . ( 1551) 
we get 

Al2 ~ 6 (e^ + iee-*) ' 11 ~ p(e^ + i«r*) [W) 
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The one-soliton solution is u(x,t) = iA^/A^, and taking into account the property 

2i cos 2 Pi 



A^ = A 1 ^*, we have 



u(x, t) = p 



(61) 



e sinh(z + iPi) + i 

The case e = —1(1) corresponds to bright (dark) soliton. The parameters k and v 
in (1591) are the soliton inverse width and the soliton velocity, respectively. In fact, 
there is only one parameter Pi characterizing the soliton, and it is usually called a one- 
parametric soliton [5]. Amplitudes (with respect to the background) of the bright and 
dark solitons are Aj, = 2psin/?i and Ad = p— p\l — 2sin/3i|, respectively. Dependences of 
the amplitudes on the parameter P are shown in figure [T](a). It is interesting to note that 
the dark soliton amplitude is a nonmonotonic function of P and the maximum occurs 
at P cr = 7r/6. The dark soliton profiles \u(z)\ for different P and p — 1 are presented in 
figure QJb). The dark soliton with P = P cr (the curve 1 in figure [Qb)) may be called 
"black" soliton: the intensity in the center of the soliton falls to zero. The corresponding 
one-soliton Jost solutions can easily be obtained from equations (jUjl . (jl6l) and (1561 

M*,C,t) = jj ppsinA + (C 2 - p 2 )A { £]h((,t), (62) 

^i(^Ct) = c( e c2 _ c ! 2) [VsinA + (C 2 - p 2 )A^]h-\C,t), (63) 



= (c2 _ c , 2 \ [2psinA + (C 2 - ^MSJWC,*), (64) 



fcfoC,*) = * ^ [2»p 2 S in/?i + (C 2 - p 2 )4 2 ) ]^ 1 (C,t). (65) 

The remaining Jost solutions can be found from the symmetry properties (I2T1) . (|22|) and 
(El 



^2 = $L> ^2 = ~$L, ¥?2 = ^J, <^2 = (66) 

Next we write out solutions for two more cases: the case when a(() has two simple zeros 
in the first quadrant on the p circle (i. e. M = 2, iV = 0) so that 

Ci = pexp(iPi), (2 = pexp(z/5 2 ), (67) 

and the case when a(() has one simple zero in the first quadrant outside the p circle (i. 
e. M = 0, N = 1) 

Ci = pexp( 7l + z/3 1 ), 7l >0. (68) 

The case ( IBTl) corresponds to two-soliton solution for the one-parametric solitons, while 
the case (|68l) corresponds to two-parametric one-soliton solution. In both cases we need 
to solve a system of four linear algebraic equations. Under this, the corresponding 
minors and determinants can be factorized and some parts of them are cancelled so that 
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the resulting expressions for u can significantly be simplified. The solutions are of the 
form 

BD , , 

« = Pd*> ( 69 ) 

where for the case (E 



B = \- te ie - 3i ^ - ie 2 e-^~ z * - e& ^f 1 ^] e -HPi+M e - Zl - z ^ (7Q) 

sm (Pi + fo) 

D = 1 — iexe^-' 1 - ie 2 e^~ z * - e ^l^ 1 ~ ^\ e^ + ^ e~ z ^ z \ (71) 

sm (/3i + (3 2 ) 

where z = k 0i j(x — Vjt — x j) (j — 1, 2) with 

fcoj = p 2 sin(2/^), Vj = 2p 2 - p 2 cos(2^), (72) 

and, as before, 6j = —1(1) corresponds to bright (dark) soliton. Equations (1691 . (1701 
and (j7ip describe collisions between bright/dark and bright/dark solitons. 
For the case (1681) we get 

£> = sinh 27! cosh(z + 271 + 3ifii — In sinh 271) + sin 2/3i sinh(37i — if), (73) 

D = sinh 271 cosh(,2 + 271 — if3\ — In sinh 271) — sin2/3i sinh(7i + iip), (74) 

with 

z = k (x — vt — x ), ip = p(x — wt) + <f , (75) 

k = p 2 cosh 27! sin 2/3 1; p = p 2 sinh 271 cos 2/?i, (76) 

u = 2p -p cos2/?i — — — , w = 2p -p cosh 271 — . (77) 

cosh 27! cos 2pi 

The two-parametric soliton given by (1691) . (1731 and (1741) with the parameters 71 and (3\ 
is actually a breather (oscillating soliton) with period 

rp (78) 

p 2 tanh(27 1 )[cosh 2 (2 7l ) + cos 2 ^)] ' 1 ' 

and with velocity v given by ( 1771) . If 71 — > and <po 7^ nit (n is an integer), we 
have T — > 00 and the breather reduces to the one-parametric soliton (bright or dark, 
depending on <p ) given by (foTi) . The found soliton solutions perfectly coincide with 
those obtained in [T5J ITT] . 

4. Integrals of motion 

Being completely integrable, the DNLSE with NVBC has an infinite set of integrals of 
motion. Eliminating (p 2 from ([6]), and substituting 

ifi = exp{— id — ik(()x + ir]~(x) + q(x, ()}, (79) 

into the resulting equation for ipi, we get Riccati equation for the function / = 
i{p 2 -\u\ 2 )/2 + d x q 

d x f + (/ - ik) 2 - ^(/ - ik) + A 2 (\u\ 2 + A 2 ] = 0, (80) 
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where v! = d x u. Representing 

«*.o=7i:^. pi) 

n=0 S 

and substituting ( 15T1) into ( 1501) . one can successively determine the coefficients f n (x). 
The first few of them are 

/o = ^(M 2 -P 2 ), (82) 

A = -m(9^*-i(M 4 -p 4 ), (83) 

f 2 = 2i(\u\ 2 - p 2 )ud x u* + i\u\ 2 u*d x u - 2ud 2 x u* + \u\\\u\ 2 - p 2 ). (84) 
From equations ({TBI) and (fT5l) we have 

y?i -> e ^- ifcx , as x -> -oo (85) 

V?i -> ae^ fcx ' - llhe ie+ikx , as a; -> oo. (86) 

It then follows from r?~(— oo) = and from equation ([79]) that g(— oo,£) = 0. Since 
?? - (oo) = 77, from equations (1791) and (1561) one can find that g(oo, Q = lna(C) + 2i0 — ir] 
as x —>■ 00 and |£| — > 00. On the other hand, from the definition of the function f(x, () 
we have 



9(00,0=/ f(x,C)dx-irf. (87) 
./ —00 

Thus, taking into account (IHTj) . one obtains 

00 r 

lno(C) = -2ie-< 53^, (88) 



n=0 

where 



/n= / /n(«)dE (89) 
J — oo 

are integrals of motion. As usual, expanding (1371) in power series with respect to l/£ 
and using (1551) . one can explicitly express the integrals of motion in terms of discrete 
(solitonic) and continuous scattering data. In particular, for I we get equation 
and for Ii we have 



A = ;£[/> 4 (i-^) + (cr-d+;f;(c 

i=l L v ^ ^ 7 J fc=l 



M 



, p\n{l-\b{p)\ 2 )dp. (90) 
5. Perturbation theory 

In the presence of perturbations the DNLSE can be written as 

id t u + d 2 u + id x (\u\ 2 u) = p[u, «*], (91) 
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where the perturbation is represented by the term p[u, u*]. Equation (I9~TI) can be cast 
in the matrix form 

d t U - d x V + [U,V] + P = 0, (92) 

where 

P=( ° iXp ). (93) 

I -iXp* J K ' 

Then, evolution equation for the monodromy matrix S can be obtained in a way similar 
to that described in [21]. As a result, we have 

POO 

d t S(t,C)+M(C)[<T 3 ,S(t,0] = - / (M + )- 1 (x,t,C)PM-(x,t,C)dx. (94) 



The equations of motion for the coefficients a(t, £) and b(t, () are contained in Eq. ( 194"1) . 
Taking into account that detM 1 * 1 = 1 + p 2 /( 2 and equation (TT2"]) . we have 

da i((p 2 - C 2 ) f 00 

(pip2(p2+p*ip 1 <pi)dx, (95) 



dt + C 2 ) 

^ - 2zO(C)6 = - -JpT^f} f (p&¥* + dx. (96) 

The expression defining the zeros Q(t) of a(t, Q is a(t,Q{t)) — 0. Differentiating with 
respect to t gives 



dC- 

d t a(t,C j (t)) + -j£a' j = 0, (97) 
where a'- = da/d(\^ = ^. Using fl95l) and fl97l) . one therefore finds 



m ~ (p 2 + c 2 K 



(P^2,jf2,j + p*^i,jipi,j) dx, (98) 



or, taking into account (l2"7]l . 



where ip2,j, <f2,j, i>i,j, and <pij are the corresponding Jost solutions evaluated at £ = Q. 
Evolution equation for bj can be obtained in a way similar to that described in Ref. 
[2Tj . As a result, one obtaines 

+P>i^7 ((pi -b^ t ) \dx', (100) 



where, after differentiating, the integrand is evaluated at ( = Q. Equations fl95l) . fl96|) . 
fl98|) and (11001) describe the evolution of the scattering data. 

If p[u, u*] is a small perturbation, one can substitute the unperturbed iV-soliton Jost 
solutions if), ip, and <p into the right-hand side of (|95|) . (|96|) . (I9~8l) and fllOOl) . This yields 
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Figure 2. Dependence of the soliton parameter (3(t) on time for different initial j3(0) 
and for (a) bright and (b) dark solitons. 



evolution equations for the scattering data in the lowest approximation of perturbation 
theory. This procedure can be iterated to yield higher orders of perturbation theory. 
The appearing hierarchy of equations is applied to an arbitrary number of solitons and, 
in particular, describes nontrivial many-soliton effects in the presence of perturbations. 
In this paper we restrict ourselves to the case of one-parametric one-soliton solutions 
with (i = pexp(i/3i) and substitute unperturbed one-soliton Jost solutions ( |63i) -(l66l) 
into the right-hand side of (195]) . (I96I) . (1981) and fllOOl) . The resulting equations are the 
desired set describing the evolution of the scattering data (both solitonic and continuous) 
in the presence of perturbations. Under this, equations fl98|) and fllOOp corresponds to 
the so called adiabatic approximation, when an unperturbed instantaneous shape of 
one soliton with slowly varying parameters j3i and bi is assumed, while equations (1951) 
and (}96l) account for radiative effects. Making use of the relation between the soliton 
solution (EH) and the corresponding squared Jost solution evaluated at f3\, the adiabatic 
equation for f3\ can be simplified to 

^ = l - J°°yu s - K) dx, (ioi) 

where u s is the one-parametric soliton solution (IMj) . Note, that this equation can also 
be obtained with the aid of the integral of motion I . 
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6. Application 

As an example of using of the present perturbation theory, we consider the case when 
the perturbation term p in (19T1) has the diffusive form 

V-iD^. (102) 

This form of dissipative perturbation occurs for Alfven solitons in a plasma when finite 
electric conductivity (and/or ion viscosity) of the plasma is taken into account [6l |2T]. 
The conditions (in terms of the plasma parameters) under which the diffusive term 
( 11021) can be considered as a small perturbation are given in [61 [21] . We consider the 
action of perturbation on the one-parametric soliton flBTl) in the adiabatic approximation. 
According to this approximation, the parameter (3 of the soliton ( 16TI) is considered as 
slowly varying in t but with the unchanged functional shape. Then, substituting (11021) 
into ( 11011) and calculating integrals with u s given by ( J6TT) . one can obtain 
36 

= -4Dp 2 sin (3[e sin (3(cos 2 (3 -3)(n -2e(3) + 2 cos (3(3 -2 cos 2 (3)} (103) 

Numerically found solutions of (11031) for different initial values of (3 are shown in 
figures [2](a) and [2](b) for bright (e = —1) and dark (e = 1) solitons respectively. 
For sufficiently small initial j3(0) <C 1, from ( 11031) one can get a simple estimate 
P(t) = (3(0) exp(— 8Dp 2 t) both for bright and dark solitons, so that their amplitudes 
and velocities decrease with time. The situation, however, changes dramatically when 
(3(0) is not too small. Under this, the behaviour of bright and dark solitons is essentially 
different. First of all, as one can see in figure [21 dark solitons turn out to be much more 
robust. Next, if the initial (3 exceeds the critical value (3 cr = n/6 (see figure CP, then 
the amplitude A d (with respect to the background) of the dark soliton first increases 
with time, reaches a maximum for (3 cr = n/6, where Ad = p, and finally decreases. 

7. Conclusion 

We have presented a simple approach for finding iV-soliton solution and the 
corresponding Jost solutions of the DNLSE with NVBC It is important that the 
exact solutions can be obtained without explicit determining of the phase factor. The 
found one- and two-soliton solutions perfectly coincide with those obtained in [151 fTTj . 
but, unlike [I5j [T7], our method allows to get solutions describing collisions between 
breathers, as well as collisions between pure bright/dark solitons and breathers. 

We have also developed a perturbation theory based on the 1ST for perturbed 
DNLSE solitons. This approach fully uses the natural separation of the discrete and 
continuous degrees of freedom of the unperturbed DNLSE with NVBC. We have derived 
evolution equations for the scattering data (both solitonic and continuous) in the 
presence of perturbations. As an application of the developed theory, we considered (in 
the adiabatic approximation) the action of the diffusive-type perturbation on a single 
bright/dark soliton. 
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